On primes p for which d divides OTdp{g) 



^ ■ Pieter Moree 



■ Abstract 



o 

s 



Let Ng{d) be the set of primes p such that the order of g modulo p, ordp{g), 
is divisible by a prescribed integer d. Wiertelak showed that this set has 
a natural density, 5g{d), with 5g{d) £ Q>o- Let Ng{d){x) be the number 
of primes p < x that are in Ng{d). A simple identity for Ng{d){x) is 
established. It is used to derive a more compact expression for 5g{d) than 
known hitherto. 



^ 1 Introduction 

(N 

I Let 5^ be a rational number such that g ^ {—1,0, 1} (this assumption on g will 

^ I be maintained throughout this note). Let Ng{d) denote the set of primes p such 

'nI" ■ that the order of g{mod p) is divisible by d (throughout the letter p will also be 

used to indicate primes). Let Ng{d){x) denote the number of primes in Ng{d) 
not exceeding x. The quantity Ng{d){x) (and some variations of it) has been the 
subject of various pubUcations [1, 3, 4, 7, 9, 11-19]. Hasse showed that Ng{d) 
has a Dirichlet density in case d is an odd prime ||3j, respectively (i = 2 jlj. The 
^ ■ latter case is of additional interest since Ng{2) is the set of prime divisors of the 

^ . sequence {g^ + Ij^li. (One says that an integer divides a sequence if it divides 

at least one term of the sequence.) Wiertelak established that Ng{d) has a 
natural density 6g{d) (around the same time Odoni P did so in the case d is a. 
prime). In a later paper Wiertelak [15] proved, using sophisticated analytic tools, 
the following result (with Li(a;) the logarithmic integral and with uj{d) = J2p\d 
which gives the best known error term to date. 

Theorem 1 [T3]. We have 

Ng{d){x) = 5gid)Uix) + Od,g f --^(loglogx)-('^)+i 

ylog X 

Wiertelak also gave a formula for 6g{d) which shows that this is always a positive 
rational number. A simpler formula for Sg{d) (in case g > 0) has only recently 
been given by Pappalardi TD]. With some effort Pappalardi's and Wiertelak's 
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expressions can be shown to be equivalent. 

In this note a simple identity for Ng{d){x) will be established (given in Propo- 
sition Q). From this it is then inferred that Ng{d) has a natural density Sg{d) that 
is given by (jU, which seems to be the simplest expression involving field degrees 
known for Sg{d). This expression is then readily evaluated. 

In order to state Theorem |21 some notation is needed. Write g = ig^, where 
go is positive and not an exact power of a rational and h as large as possible. 
Let D{gQ) denote the discriminant of the field Q(y^). The greatest common 
divisor of a and h respectively the lowest common multiple of a and h will be 
denoted by respectively [a, 6]. Given an integer d, we denote by the 

supernatural number (sometimes called Steinitz number), lipid Note that 

{v.d°^) = Y[^^,p'^^^^\ 

Definition. Let d be even and let eg{d) be defined as in Table 1 with 7 = 
max{0, U2{D{go)/dh)}. 

Table 1: eg{d) 



9\l 


7 = 


7 = 1 


7 = 2 


g>0 


-1/2 


1/4 


1/16 


g<0 


1/4 


-1/2 


1/16 



Note that 7 < 2. Also note that eg{d) 
Theorem 2 We have 



1/2)2" if > 0. 



Sg{d) 



ei 



n 



d(h,d°°) ^^p- 



with 



( 1 



ei 



l + 3(l-sgn((7))(2 



l)/4 



if2\d; 
if 2\ \d and D{gQ 



Ad; 



l + 3(l-sgn(^7))(2'^2W-l)/4 + e,(rf) if 2\\d and D{go)\4d; 



1 



^l + eigi{d) 
In particular, if g > 0, then 

1 

and if h is odd, then 



ei 



ei 



-1/2) 



2max{0,^2{D{9)/dh)} 



tf4\d, D{go)^4d; 
tfA\d, D{go)\Ad. 



if2\d and D{gQ)\Ad; 
otherwise, 



if2\d and D{g)\4d; 
otherwise. 



Using Proposition n of Section 2 it is also very easy to infer the following result, 
valid under the assumption of the Generalized Riemann Hypothesis (GRH). 

Theorem 3 Under GRH we have 

Ng{d){x) = 6g{d)U{x) + Od,g(v^log"('')+^ x), 
where the implied constant depends at most on d and g. 

In Tables 2 and 3 (Section 6) a numerical demonstration of Theorem El is given. 
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2 The key identity 



Let ttl{x) denote the number of unramified primes p < x that spht completely 
in the number field L. For integers r\s let Kg ^ = QiCs, g^^"^)- 

The starting point of the proof of Theorem |21 is the following proposition. By 
rp{g) the residual index of g modulo p is denoted (we have rg{p) = [¥p : {g)]). 
Note that oTdp{g)rp{g) = p — 1. 

Proposition 1 We have Ng{d){x) = J2v\d-- EaM/^(«)^ifd.,«. 



X 



Proof. Let us consider the quantity Xlaid /^('^)^i^d„ ^ prime p counted 

by this quantity satisfies p < x, p = l(mod dv) and rp{g) = vw for some in- 
teger w. Write w = W1W2, with wi = {w,d). Then the contribution of p to 
J2a\df^i^)^K,,,c,S^) is Ea|«;i/^(")- We couclude that 

a\d 

< X : p = l(mod dv), v\rp{g) and {— , d) = 1}. (1) 

It suffices to show that 

Ng{d){x) = ^ #{j9 < X : J9 = l(mod dv), v\rp{g) and (^^, d) = 1}. 

v\d°° 

Let p be a prime counted on the right hand side. Note that it is counted only once, 
namely for v = {rp{g),d°°). From oTdp{g)rp{g) = p — 1 it is then inferred that 
d\oTdp{g). Hence every prime counted on the right hand side is counted on the 
left hand side as well. Next consider a prime p counted by Ng{d){x). It satisfies 
p = l(mod d). Note there is a (unique) integer v such that f p = l(mod dv) 
and {rp{g)/v, d) = 1. Thus p is also counted on the right hand side. □ 

Remark 1. From (Q) and Chebotarev's density theorem it follows that 

< y ^^"^ < I . (2) 

" ^ [Kdv,av ■■ Q] - [Kdv,v : Q] ^ ' 

3 Analytic consequences 

Using Proposition Q it is rather straightforward to establish that Ng{d) has a 
natural density 6g{d). 

Lemma 1 Write g = gi/g2 with gi and g2 integers. Then 

'(log logx)'^^'*)' 
log^/® X 



Ng{d){x) = [ Sg{d) + Od,g { ' t ) ) Li(a:), (3) 



where the implied constant depends at most on d and g and 



v\d°° a\d 
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Corollary 1 The set Ng{d) has a natural density 5g{d). 

The proof of Lemma Q makes use of the following consequence of the Brun- 
Titchmarsh inequality. 



Lemma 2 Let 7r(x; /, k) = ^p<^^ ^^^^^^^ 1. Then 

} 7r(x; dv, 1) =0d[ 



ui(d) 



uniformly for 3 < z < ^Jx. 

Proof On noting that Mii{x) := #{f < x : v\d°°} < (logx)'^*^'^Ylos2, it straight- 
forwardly follows that 



^ V 

>2 



dMd{z) (log^)-W 



By the Brun-Titchmarsh inequality we have 7i{x;w,l) <^ x/{ip{w)\og{x/w)), 
where the implied constant is absolute and w < x. Thus 



7r(x; , 1) <C 



5^ - <d 



X (log z)'^^'^^ 



z<v, dv<x^/^ 



ip{d) logx ^ t; logx z 



(5) 



1) > z 



Using the trivial estimate 7r(x; d, 1) < x/d we see that 



V vr(x;rft;,l)< V f < V - «, x^/=^(logxr(^). 



On combining (jSj) and © the proof is readily completed. 



(6) 



□ 



Proof of Lemma^ From [lOj Lemma 2.1] we recall that there exist absolute 
constants A and B such that if f < BilogxY^^ /d, then 



dv .av 



Lifx) 



+ Oo(xe~ 



dv,av 



(7) 



Let y = B(\ogxY^^ / d. From the proof of Proposition Q we see that 



Ng{d){x) = 'Yfi{a)7iK^^ ^^{x) +0 7i{x;dv, 1) 



a\d 

v<y 



v>y 
\v\d°° 



h + Oih), 



J 



say. By LemmaElwe obtain that I2 = 0(x(loglogx)'^'^'^) log x). Now, by (jZj), 
we obtain 



"1^°° Q\d 



dv,av 



X 



Denote the latter double sum by I^- Keeping in mind Remark 1 we obtain 



h = S,{d) + 



EE 



\ 



1-1^°° a\d 

v>y 



dv,av 



Using (j21) and Lemma El it follows that 

/ 



EE 

a\d 

v>y 



O 



dv,av ■ 



E — 



v>y 



/i(log|/)-('^) ^(logy)-W 

^d[ ) — 'Jd,g 



y 



y 



and hence 



h = Sg{d) + Od,g 



(logy) 



ui(d) 



y 



The result follows on collecting the various estimates. 



□ 



4 The evaluation of the density Sg{d) 

A crucial ingredient in the evaluation of Sg{d) is the following lemma. 

Lemma 3 [6J. Write g = ^gQ, where g^ is positive and not an exact power of 
a rational. Let D{gQ) denote the discriminant of the field Q.i^/gt))- Put m = 
D{go)/2 ifP2{.h) = andD{go) = 4(mod 8) or i^2ih) = 1 and D (go) = O(mod 8), 
and m = [2'^^^^^^'^, D^go)] otherwise. Put 

( m if g < and r is odd; 

nr= <^ p^2{hr)+i^ otherwise. 

We have 

where, for g > or g < and r even we have 

^ ^ [1 ifnr\kr, 
and for g < and r odd we have 

{2 ifurlkr; 
\ zf2\k and2''^^''^+^ ] k; 
1 otherwise. 

Remark. Note that if h is odd, then Ur = [2'^^^'^^+^^ Did)]- Note that Uj. = n^^ir)- 

The 'generic' degree of [Kdv,av '■ Q] equals ip{dv)av/ {av, h) and on substituting 
this value in (jH) we obtain the quantity 5*1 which is evaluated in the following 
lemma. 
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Lemma 4 We have 

(p{dv)av 



aid 

where 



p\d 

Proof. Since for ^10?°° we have ip{dv) = v(f{d), we can write 



v\d°° a\d ^ ' v\d'^ Q.\d 



(fid) ^ ^ av^ ip(d) ^ ^ a(v,h) 

' v\d'=^ o\d ' v\d'=^ a\d y ' I 

The expression in the inner sum is multiphcative in a and hence 



^ if 

otherwise. 



^ li{a)[av,K) ^ TT A _ 
^ a(v,h) AAl p(v,h)l 

a\d ^ ' ^ p\d ^ ^\ ^ J / 

On noting that {v, h)/v^ is multiphcative in v and that for k > ^^{h) 

oo 

E 



{p\ h) _ pMh)+2-2k 

r=k 



p2r — 1 ' 

one concludes that 



This completes the proof. □ 

Remark. Note that the condition (/i, 0?°°) is equivalent with fpiv) > Pp{h) 

for all primes p dividing d. 

By a minor modification of the proof of the latter result we infer: 
Lemma 5 Let k > be an integer. Then 

v\d°° a\d ' 

The next lemma gives an evaluation of yet another variant of S\. 
Lemma 6 Let D he a fundamental discrimant. Then 



"1^°° o.\d ^ ' I U otherwise, 

[2'"2(''d/a)+l 

where 7 = max{0, 1/2(0 /dh)}. 



4:~^S{d, h) if2\d, D\4:d and 7 > 1; 



Proof. The integer ^2'^'^^'^'^^"^^^, D] is even and is required to divide d°°, hence 
5*3 (D) = if d is odd. Assume that d is even. If D has an odd prime divisor not 
divinding d, then D \ d°° and hence S^lD) = 0. On noting that V2{D) < 1^2(4^) 
and that the odd part of D is squarefree, it follows that if 5*3(1^) 7^ 0, then 
D\4:d. So assume that 2\d and D\4:d. Note that the condition [2"^'^'''^/°'^+'^ , D]\dv 
is equivalent with z/2(f ) > z/2(/?.) + max{l, h'2{D/dh)} for the a that are odd, and 
M'^) > Mh) + 7 for the even a. Thus if 7 > 1 the condition [2'^2(/^rf/a)+i^ jji^^^y 
is equivalent with i'2{v) > i'2{h) +7 and then, by Lemma El S^{D) = 5*2(7) = 
A-^Sid,h). If 7 = then 



S,{D) = S2iO)- E 



ip(dv)av 

ti|dO° a\d ~ ^ ' 



By Lemma El it follows that 5*2(0) = Sid^ h). A variation of Lemma 01 yields that 
the latter double sum equals 3S{d, h)/2. □ 

Remark. Put 

e^iD) = I (-1/2)2'"^''"'"'^'°'*'^ if 2\d and D\4d; 
1 otherwise. 

Note that Lemma IHl can be rephrased as stating that if D is a fundamental dis- 
criminant, then S^i^D) = e2{D)S{d, h). 

Let g > 0. It turns out that ordp(5') is very closely related to ordp(— (7) and this 
can be used to express N_g{d){x) in terms of Ng{*){x). From this S^g{d) is then 
easily evaluated, once one has evaluated Sg{d). 

Lemma 7 Let g > 0. Then 

/V (.^(.^fNgil)ix)+Ngi2d)ix)~Ngid)ix) + 0il) ^/ ^ 2 ( mo d 4 ) 
\n^^cI){x) + 0{1) otherwise. 



In particular, 

5.g{d) 



6g{f, + 5g{2d) ~ 5g{d) if d = 2{mod A) 



Sg{d) otherwise. 

The proof of this lemma is a consequence of Corollary ^ and the following obser- 
vation. 

Lemma 8 Let p be odd and g ^ be a rational number. Suppose that Vpig) = 0. 
Then 

( 2ordp(^) z/2tordp(^7); 
OTdp{-g) = I oTdp{g)/2 if ordp{g) = 2(mod 4); 
[oTdp{g) if 4:\oTdp{g) . 

Proof. Left to the reader. □ 

Remark. It is of course also possible to evaluate Sg{d) for negative g using the ex- 
pression and Lemma El however, this turns out to be rather more cumbersome 
than proceeding as above. 
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5 The proofs of Theorems |2l and 131 

Proof of Theorem^ By Lemma [T] it suffices to show that 

v\d°° a\d 

If g > 0, then it follows by LemmaElthat Sg{d) = 5*1 + 6*3 ((70)) and by Lemmas 
13] and ini (with D = D{gQ)), the claimed evaluation then results in this case. If h 
is odd, then similarly, 6g{d) = Si + 5*3 (cf. the remark following Lemma 
IHl) and, again by Lemma 0] and the claimed evaluation then is deduced in this 
case. If g < 0, the result follows after some computation on invoking Lemma 
and the result for g > 0. □ 

Proof of Theorem^ Recall that 71l{x) denotes the number of unramified primes 
p < X that split completely in the number field L. Under GRH it is known, cf. 
0, that 



where di denotes the absolute discriminant of L. From this it follows on using 
the estimate log \dK^^_^ < dv{\og{dv) + log \gig2\) from [H] that, uniformly in v, 

[J^dv,av ■ m 

where a is an arbitrary divisor of d. On noting that in Proposition ^ we can 
restrict to those integers v satisfying dv < x and hence the number of non- 
zero terms in Proposition ^ is bounded above by 2'^^'^-*(logx)'^'''^'', the result easily 
follows. □ 



6 Some examples 

In this section we provide some numerical demonstration of our results. 

The numbers in the column 'experimental' arose on counting how many primes 
P < Pio» = 2038074743 with I'p^g) = 0, satisfy d\oTdp{g). 



Table 2: The case > 



9 


^0 


h 


Digo) 


d 


ei 


5M 


numerical 


experimental 


2 


2 


1 


8 


2 


17/16 


17/24 


0.70833333 ■ ■ ■ 


0.70831919 


2 


2 


1 


8 


4 


5/4 


5/12 


0.41666666- ■■ 


0.41667021 


2 


2 


1 


8 


8 


1/2 


1/12 


0.08333333 ■ ■ ■ 


0.08333144 


3 


3 


1 


12 


11 


1 


11/120 


0.09166666- ■■ 


0.09165950 


3 


3 


1 


12 


12 


1/2 


1/16 


0.06250000 ■ ■ ■ 


0.06249098 


4 


2 


2 


8 


5 


1 


5/24 


0.20833333 ■ ■ ■ 


0.20833328 


4 


2 


2 


8 


6 


5/4 


5/32 


0.15625000- •■ 


0.15625824 
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Table 3: The case < 



9 


9q 


h 


^(^o) 


d 


ei 


6g{d) 


numerical 


experimental 


-2 


3 


1 


8 


2 


17/16 


17/24 


0.70833333 ■ ■ ■ 


0.70835101 


-2 


2 


1 


8 


4 


5/4 


5/12 


0.41666666- ■ ■ 


0.41667021 


-2 


2 


1 


8 


6 


17/16 


17/64 


0.26562500- ■■ 


0.26562628 


-3 


3 


1 


12 


5 


1 


5/24 


0.20833333 ■ ■ • 


0.20834107 


-3 


3 


1 


12 


12 


1/2 


1/16 


0.06250000- ■■ 


0.06249098 


-4 


2 


2 


8 


2 


2 


2/3 


0.66666666 ■ ■ ■ 


0.66666122 


-4 


2 


2 


8 


4 


1/2 


1/8 


0.08333333 ■ ■ ■ 


0.08333144 


-9 


3 


2 


12 


2 


5/2 


5/6 


0.83333333 ■ ■ ■ 


0.83333215 


-9 


3 


2 


12 


6 


11/4 


11/32 


0.34375000 ■ ■ ■ 


0.34375638 



Acknowledgement. I like to thank Francesco Pappalardi for sending me his 
paper [10- Theorem 1.3 in that paper made me realize that a relatively simple 
formula for 5g{d) exists. The data in the tables are produced by a C^^ program 
kindly written by Yves Gallot. 
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